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We analyzed the 100-yr light curves of Galactic high-mass X-ray binaries using the Harvard pho- 
tographic plate collection, made accessible through the DASCH project (Digital Access to a Sky 
Century at Harvard). As scanning is still in progress, we focus on the four objects that are cur- 
rently well covered: the supergiant X-ray binary Cyg X-l (VI 357 Cyg), and the Be X-ray bina- 
ries 1H 1936+541 (BD+53 2262), RX J1744.7-2713 (HD 161103), and RX J2030.5+4751 (SAO 



o 

49725), the latter two objects being similar to y Cas. The star associated with Cyg X-l does 
not show evidence for variability with an amplitude higher than 0.3 magnitude over a hundred 
1-1 years. We found significant variability of one magnitude with timescales of more than 10 years 

^-H for SAO 49725, as well as a possible period of ~500-600 days and an amplitude of 0.05 mag- 

nitude that might be the orbital, or super-orbital period of the system. The data is insufficient to 
conclude for HD 161103 but suggests a similar long-term variability. We thus observe an addi- 
tional characteristic of y Cas-like objects: their long-term variability. This variability seems to 
be due to the slow evolution of a decretion disk around the Be star, but may be triggered by the 
presence of a compact object in the system, possibly a white dwarf. This characteristic could be 
used to identify further similar objects otherwise difficult to detect. 
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1. Introduction 

High-mass X-ray binaries (HMXBs) are composed of a compact object (a neutron star or 
a black hole) that accretes material from a massive early-type O-B star. Conventionally, they are 
subdivided into two groups, the supergiant X-ray binaries (sgXBs) and Be/X-ray binaries (BeXBs). 

In sgXBs, the massive star (type I/II star) ejects a slow, dense wind and the compact object 
directly accretes the stellar wind, leading to a persistent X-ray emission (Lx ~ 10 35-36 erg s _1 ). 
They exhibit rare Type II outbursts and no Type I outbursts. Some sgXBs are Roche lobe over- 
flow systems causing a high X-ray luminosity (L x ~ 10 38 erg s _1 ) during outbursts. In particular, 
Cyg X-l is the only sgXB with Roche lobe overflow and stellar wind accretion hosting a confirmed 
black hole. It is generally found in the low state. 

A BeXB consists of a neutron star orbiting a Be star in a wide (period of months) and eccentric 
orbit. The Be star is surrounded by a disc of relatively cool material, presumably ejected from the 
star due to causes unknown, but generally believed to be associated with fast rotation, magnetic 
fields and/or non-radial pulsations [1]. In BeXB systems, the neutron star can penetrate this equa- 
torial disc during each periastron passage, giving rise to periodic outbursts (Type I) that can be 
observed over a wide range of wavelengths (optical, X-ray, ...). Be stars are variable and may lose 
their disk (and hence emission) completely and be classified as normal B-type stars. Therefore, 
census studies will provide a lower limit to the fraction of Be stars within a sample. A scenario 
for B stars to become Be stars proposes the rapid rotation of Be stars as consequence of binary 
evolution spin-up [2], Hence all Be stars would have evolved companions. The Be star is hardly 
affected in a BeXB: while the disks' radii have an effective upper size, probably determined by the 
binary parameters [3], they are otherwise not significantly different from single Be star disks. A 
recent review of basic properties of BeXB can be found in [4]. 

A growing number of early Be stars discovered in X-ray surveys with low-luminosity hard 
X-ray emission (Lx ~ io 32-33 erg s _1 ) show a very likely thermal emission of high temperature 
(kT ~ 8 keV). A strong iron line, including a fluorescence component, was detected for some 
systems. They also share very similar optical properties such as spectral types (B0.5III-Ve) and 
strength of the Balmer emission. This pattern of common X-ray and optical properties appears 
quite similar to that of y Cas and points to the emergence of a new class of y Cas analogs [5]. Two 
models have been put forward to explain this type of system: accretion onto a compact object (most 
likely a white dwarf) and magnetically heated material between the photosphere of the B star and 
the inner part of its disc [6]. 

The Digital Access to a Sky Century at Harvard (DASCH) is a project to digitize and analyze 
the scientific data contained in the ^530000 Harvard College Observatory (HCO) plates taken 
between the 1880s and 1990s, which is a unique resource for studying temporal variations in the 
universe on 10-100 yr timescales [7]. We have developed the astrometry and photometry pipeline 
[8, 9, 10]. We started to scan plates in several selected fields and discovered new types of variable 
stars [11, 12, 13]. 

We initiate here the study of the long-term variability of known HMXBs in the Galaxy by 
studying the first four objects covered by the current DASCH scanned images. DASCH will enable 
a large study of the long-term variability of Galactic HMXBs and the brightest HMXBs in the 
Magellanic Clouds. 
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X-ray source name Star name 



Spectral type Bmag 



RA Dec 



RX J1744.7-2713 HD 161103 

1H 1936+541 BD+53 2262 

Cyg X-l / 4U 1956+3 V1357 Cyg 

RX J2030.5+475 1 SAO 49725 



B0.5III-Ve 9.13 



Be 10.38 
09.7Iab 9.62 
B0.5III-Ve 9.57 



17:44:45.70 -27:13:44.0 

19:32:52.30 +53:52:45.0 

19:58:21.70 +35:12:06.0 

20:30:30.80 +47:51:51.0 



Table 1: HMXB candidates with DASCH data. 



2. Data 

2.1 Search for HMXBs in DASCH data 

The DASCH database now contains measurements from 38 000 plates (November 2012), 
7% of the total collection, covering mainly five selected fields in the sky. As an initial study 
of the HMXB population, we used the catalog of HMXBs from [14] and searched DASCH for 
matches. This catalog contains 114 HMXBs, however, some sources are only tentatively identified 
as HMXBs on the basis of their X-ray properties similar to the known HMXBs. About 60% of the 
catalogued HMXB candidates are known or suspected BeXB, while 32% are sgXB. 

We found that four objects have DASCH light curves with more than 500 magnitude measure- 
ments over 100 years. They are listed in Table 1 with their spectral type, coordinates and current 
B magnitude (Bmag) as found in SIMBAD 1 . One of those objects is the sgXB Cyg X-l, and the 
three others are BeXB. We note that two objects (HD 161 103 and SAO 49725) are peculiar HMXB 
candidates similar to y Cas for which the binary nature remains to be proven [5]. 

2.2 Cleaning of the light curves 

We removed all points from plates flagged as showing defects, and we performed a 3d clipping 
to remove isolated outliers in order to focus on long-term variability. Points flagged as saturated are 
plotted with crosses and might be affected by calibration issues, though they may contain relevant 
information. The four light curves listed in Table 1 are shown in Figures 1, 2, 3, and 4 (left panel). 

Cyg X-l = V1357 Cyg is close to V1674 Cyg (55"). Both are variable stars with B mean 
magnitudes 9.62 and 10.6, respectively. For some plate series with low resolution (6" per pixel), 
the two stars can be blended. The brightest magnitudes observed in the uncleaned light-curve (9.0) 
matches well the expected magnitude of the two stars blended together, raising doubts on any real 
variability detected in the DASCH light curve. 

To complement the DASCH light curves with data after 2000, we downloaded the available 
data from the AAVSO 2 website for those objects. We used the AAVSO B magnitudes only, given 
that DASCH magnitude measurements are similar to B magnitudes. 

2.3 Periodicity search 

We only kept the non-saturated points and computed a Lomb-Scargle periodogram following 
the method described in [15]. We first subtracted the mean flux in the light curve and set the 
variance to one to obtain the normalized periodogram. We searched for periods between 100 and 

lr rhe SIMBAD database is operated at CDS, Strasbourg, France (http://simbad.u-strasbg.fr/simbad) 
2 American Association of Variable Star Observers (http://www.aavso.org) 



3 



DASCH 100-yr light curves of high-mass X-ray binaries 



Mathieu Servillat 




30000 40000 

Time (MJD) 



5000 10000 20000 




100 200 300 500 1000 2000 

Period (days) 



5000 10000 20000 



Figure 1: HD 161 103 light curve and periodograms for DASCH good points, all points (including saturated 
points shown as red crosses in the light curve) and for the detrended light curve. 
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Figure 2: BD+53 2262 light curve and periodograms. The AAVSO points are shown in green. 
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Figure 3: VI 357 Cyg light curve and periodograms. 
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Figure 4: SAO 49725 light curve and periodograms. The AAVSO points are shown in green. 
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20 000 days (~ 55 years) given the typical sampling and duration of DASCH light curves. The light 
curves and periodograms are shown in Figures 1, 2, 3, and 4 (left and right panels, respectively). 

We used a local regression technique (lowess) with a window of 7 000 days (about 20 years) 
to extract the long-term variability from the light curves. We obtained a detrended light curve after 
dividing the initial light curve by this smoothed light curve. We then recomputed a periodogram 
for the detrended light curve. 

We note that there is often a signal around 365 days (note also the weak signal at 180 days 
for HD 161103, Figure 1). This one year period seems to be due to the long-term variability of an 
unevenly sampled light curve. This phenomenon is due to spectral leakage of low-frequency power 
present in the light curves (see e.g. [16]). Indeed, the high power at one year disappears when the 
light curve is detrended, while a real signal would be strengthened. 

There is no significant variability in the light curves of HD 161 103, BD+53 2262 and V1357 
Cyg. No feature is observed in the detrended periodograms of those stars, indicating that no period- 
icity is observed. Cyg X-l has an orbital period of 5.6 days and showed ellipsoidal variations with 
a low amplitude of ^0.05 magnitude [17]. Here we can rule out modulations with an amplitude 
higher than ^0.3 magnitude. For HD 161103, we note that many points are flagged as saturated, 
but the light curve might also indicate large time scale variability. 

The variability of SAO 49725, almost one magnitude on a timescale of more than 10 years, 
is clearly seen in the light curve (Figure 4, left). We find a signal around 5 000 days (~15 years) 
in the periodogram which can be seen in Figure 4 (two cycles from 1915 to 1945). However, 
the periodicity is not seen in the first part of the light curve (1900 to 1915), indicating that this 
variability is not periodic. For the AAVSO light curve of SAO 49725, we find a possible signal 
around 500-600 days. Indeed, two cycles can be seen in Figure 4 (inset in left panel), with peaks 
around MJD 55 000 and 55 500, and an amplitude of ^0.05 magnitude. This kind of period is 
difficult to detect in the DASCH light curve due to the low amplitude and larger error bars. 

3. Discussion and conclusions 

For the first time, we study a sample of HMXBs over a hundred years using optical light curves 
from the DASCH project. This study reveals a new kind of variability, poorly studied for HMXBs. 
A more complete study will follow when the Galactic plane will be fully covered by DASCH. 

We do not expect to detect large variability amplitudes from sgXB, as illustrated by the case 
of Cyg X-l, because the optical light is mostly emitted by the supergiant star. However, Be stars 
can be extremely variable over a hundred years, as shown in Figure 4 for SAO 49725. 

It is interesting to note that among the stars presented in this study, the two sources that are 
found to be highly variable (clear variability for SAO 49725, to be confirmed for HD 161 103) are 
y Cas-like objects. This may be a remarkable property of those objects, as it reminds the long-term 
variability of y Cas itself, that varied from 1.6 to 3.0 magnitudes in the V-band between 1935 and 
1940, and then slowly brightened to 2.1 today. The long-term, non-periodic variability we observe 
for SAO 79425 could be due to changes in the decretion disk properties (size, density, temperature). 
We found a weak hint of a period around 500-600 days, consistent with an orbital or a super-orbital 
period of a BeXB (e.g. [18]). This would indicate a possible link between binarity and long-term 
variability. A similar kind of variability has been observed for some neutron stars BeXBs in the 
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Magellanic Clouds [18]. However, a clear proof of binarity is still missing for y Cas-like objects. It 
is still not clear if they are isolated Be stars, or if they host a compact object, probably a white dwarf 
to explain the hard X-ray spectrum with intermediate luminosities between Be stars and BeXB [5]. 

For the next few years, with a complete coverage of the Galactic plane, DASCH brings the 
prospect of detection of new HMXBs through the detection of long-term variability of B/Be stars. 
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